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Abstract 

The neutrino signal from SN1987A is analysed with respect to spin-flavour 
oscillations between electron antineutrinos, v e , and muon neutrinos, Vn, by 
means of a maximum likelihood analysis. 

Following Jegerlehner et al. best fit values for the total energy released 
in neutrinos, Et, and the temperature of the electron antineutrino, Tp e , for 
a range of mixing parameters and progenitor models are calculated. In par- 
ticular the dependence of the inferred quantities on the metallicity of the 
supernova is investigated and the uncertainties involved in using the neutrino 

signal to determine the neutrino magnetic moment are pointed out. 
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I. INTRODUCTION 



T 



At the end of their lives massive stars explode in a type II supernova event releasing 
almost their entire binding energies in form of neutrinos of all flavours. Depending on the 
model for the progenitor, the equation of state and the treatment of the neutrino transport 
the theoretical prediction for the total energy emitted in neutrinos is 

E t « 2 - 4 x 10 53 erg (1) 

and for the mean neutrino temperature is 

3-4 MeV for u e , 

5-6 MeV for v e , (2) 
7 — 9 MeV for all other species. 

|I],|2| The inferred values from the observed neutrino signal from the explosion of SN1987A 
in the Large Magellanic Cloud (LMC) are E t « 2-3 x 10 53 erg and T Pe « 3 MeV While 
the energy, E t , is in good agreement with supernova theory, the temperature of the electron 
anti-neutrinos lies below the predicted range. Neutrino oscillations could significantly alter 
the energy spectra as observed by neutrino detectors on Earth and have implications for 
the interpretation of the observed neutrino signal. In fact an admixture of the hotter v^- 
spectrum through oscillations would shift the inferred value for Tp e towards lower values 
thus making it even more difficult to reconcile theory with observation. 

This analysis follows closely the method by Jegerlehner et al. ||, who investigated the 
implications of flavour oscillations [p e <-> v^) on the interpretation of the supernova neutrino 
signal using a maximum likelihood method. The impact of the MSW-effect on the neutrino 
spectra has also been studied in Ref. [0-0. For a 'normal' mass hierarchy (i.e. m Ve < m u ) 
flavour oscillations between antineutrinos are effective only for a very limited choice of mixing 
parameters. However spin-flavour (SF) oscillations, which arise from the interaction of the 
(hypothetical) magnetic moment of the neutrino with an external magnetic field, could have 
a serious effect for a wide range of reasonable parameters. 



In this paper I will investigate the implications of SF-conversion between right-handed 
v e s and left-handed v^s (or v T s) on the interpretation of the neutrino signal from supernova 
SN1987A. The distortion of the electron anti-neutrino spectrum is of special interest since 
almost all the neutrinos detected from SN1987A were v e s. I will also consider the possibility 
of deriving an upper bound on the neutrino magnetic moment from the observed neutrino 
spectra. 

The conversion probabilities (i.e. the probability that a neutrino emitted from the neu- 
trinosphere has been converted into a different species by the time it leaves the supernova) 
are taken from flHj , which were calculated on the basis of the precollapse model of Woosley 
and Weaver (1995) ]ll| for a progenitor mass of 15 and 25 solar masses and solar (Z Q ) as 
well as zero (Z ) metallicity. The mass-squared difference Am 2 = m 2 ^ — m 2 e and the vacuum 
mixing angle 9 V are assumed to take the values suggested by the MSW-solution to the solar 
neutrino problem. 

In the following section I will briefly review the theory of SF oscillations. In section III I 
will discuss the various progenitor models with their associated conversion probabilities and 
in section IV I will present the numerical results which will be followed by the discussion. 



II. SPIN-FLAVOUR OSCILLATIONS 



Assuming that the neutrino is a Majorana particle with a non- vanishing transition mag- 
netic moment, an anti- neutrino of one flavour can be resonantly converted into a neutrino of 
another flavour on traversing a magnetic field B and vice versa |12|- |14]1 . These spin-flavour 
oscillations are governed by the evolution equation 



d 



dr 



V 



(3) 



\i v B 
fi„B AH 

where \x v denotes the neutrino transition magnetic moment and B the component of the 
magnetic field which is transverse to the direction of propagation of the neutrinos. The 
diagonal component of the Hamiltonian is given by 



AH = —— cos 29 v - AV (4) 
2E U 

with AV = ^/2GFp/m N (l — 2Y e ), where p denotes the density, Gp the Fermi coupling con- 
stant, 771n the nucleon mass, E u the neutrino energy and Y e the electron fraction. Resonance 
occurs when AH = and the conversion is adiabatic if 

fiB 
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dAH 



> I- (5) 



dr 

If the neutrino is a Dirac particle, the effective potential is proportional to (1 — 3Y e ), 
which is zero in the strongly neutronized core. However the oscillations are highly non- 
adiabatic in the high density core region for all but unrealistically high magnetic fields so 
that a fine tuning would be required for any oscillations to have a noticeable effect. Moreover 
SF-oscillations for Dirac neutrinos are strongly constrained by observations [|3|,|^, p~5|JT7[1 and 
subsequently it will be assumed here that the neutrino is a Majorana particle. 

It has been pointed out by Athar et al. J7| that in the isotopically neutral (i.e. number 
of protons equals number of neutrons) region above the iron core and below the hydrogen 
envelope the effective matter potential is suppressed by 3-4 orders of magnitude. Recently 



Totani and Sato |TD[ noted that the potential is even further suppressed in the O+C and 
He layers in the absence of heavy nuclei which, with their excess of neutrons, are mostly 
responsible for the deviation from isotopic neutrality and it turns out that the conversion 
probability is very sensitive to the metallicity of the supernova progenitor. 

According to the electroweak standard model the neutrino does not possess a magnetic 
moment and an indication of a non- vanishing magnetic moment would be a sign for physics 
beyond the standard model and have profound implications for particle physics. 

The current limit on the magnetic moment of the electron anti-neutrino from labora- 
tory experiments is p Ve < 2.4 x l§~ 1Q p B [050 for the Majorana neutrino. Astrophysical 



constraints on the magnetic moment from the SN1987A neutrino signal have been derived 
by Notzold [T5[ and others [IHJrTLEfl] excluding values for p v of greater than a few times 



10 pb- For a review see pl| . In the following analysis a transition magnetic moment of 



around [i v fa 10 12 /j>b will be considered. 



III. SUPERNOVA MODELS AND CONVERSION PROBABILITIES 

The value of the effective matter potential AV was calculated from the composition and 
density profiles of the precollapse model by Woosley and Weaver (1995) |IIfl , which includes 
about 200 isotopes up to 71 Ge with solar abundances. The magnitude of AV as function of 
radius is shown in Fig. l.(a) and (b) for a progenitor mass of 15 M and 25 M respectively. 
It is about 10 -5 eV in the centre and quickly decreases by about 10 orders of magnitude as 
the density decreases and Y e approaches 0.5. Then at around 0.1 R Q the effective potential 
falls discontinuously by several orders of magnitude and continues to decrease until it changes 
sign at roughly the solar radius, where the He layer ends and the hydrogen envelope starts. In 
this region Y e is very close to 0.5 (isotopically neutral) and the matter potential is suppressed 
by several orders of magnitude. In this region the deviation from Y e = 0.5 is mainly due 



to isotopes like Ne, Mg, Al, S and Ar. Totani and Sato [Kj have reduced the abundance 
of these isotopes over solar abundances thus suppressing the effective matter potential even 
further. These models, subsequently called low metallicity (Z ) models, are shown in Fig. 
1. (c) and (d). The investigation of low metallicity progenitors seems especially important 
given the low metallicity encountered in the LMC of around Z fa Z & /4 and various other 
observations that have led to the belief that the progenitor of SN1987A was a low metallicity 
blue super-giant (BSG). For a study of the progenitor of SN1987A see e.g. Langer et al. 



1989 [ 22j] . It can be seen from Fig.l that in the Zq - models the isotopically neutral region 
is broader and the effective potential in this region is suppressed by another 1-2 magnitudes 
over the solar metallicity models. 

For the purpose of calculating the conversion probabilities it can be assumed that the 
composition of the envelope remains unchanged during infall since the collapse is homologous 
and the shock wave cannot reach the isotopically neutral region during the neutrino burst 
(1-10 s). Resonance occurs in the isotopically neutral region for values of lO" 10 ^ < < 
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10 -leVi 
1U MeV 

There is no direct information on the magnitude of magnetic fields inside a collapsing 
star, but a rough estimate may be obtained from the observed field on the surface of white 
dwarfs |23f suggesting magnetic fields of up to B pa 10 10 Gauss at the surface of the iron 
core. Assuming that the radial dependence of the magnetic field above the core follows a 
dipole field, B oc r -3 , which would give a field in the isotopically neutral region of about 10 4 
G. Higher estimates may be obtained by equating the magnetic field energy to the thermal 



energy of the gas (see llQfl). In following analysis I will assume a dipole field with field 
strengths at the surface of the iron core of Bq pa 10 9 G and Bq pa 10 10 G. 

Also I assumed mixing parameters which are favoured by the MSW-interpretation of the 
solar neutrino problem: the large angle solution with Am 2 pa 10 _6 eV 2 , sin 2 26 v ~ 0.7, and 
the small angle solution with Am 2 pa 6 x 10 _6 eV 2 and sin 2 2#„ pa 10~ 2 . 

The conversion probabilities corresponding to the models discussed above have been 
calculated by Totani and Sato [llj and some typical cases are shown in Fig. 2 for Am 2 = 
10 -6 eV 2 and sin 2 2^ pa 10~ 2 . For other mixing angles or neutrino masses the energy can 
simply be rescaled with (^^7^2 cos26*„). 

These plots demonstrate that the magnitude of the conversion probability as well as 
the dependence on energy are very different for the solar and low metallicity models. For 
solar metallicity the conversion probability is generally very small for all but the most 
extreme values of B or very high values of Am 2 , which we will not consider here. For a 
field strength of Bo pa 10 9 G the conversion probability is negligible and not shown in Fig. 
2 and even for Bq pa 10 9 G the conversion probability is hardly significant. As expected 
the efficiency of the conversion improves with increasing fiB (and thus increasing 7). The 
conversion probability decreases with increasing E v as seen from the dependence of 7 on 
However for progenitors of low metallicity the conversion probability increases with 
energy and remains constant at high energies due to the precession of the magnetic moment 
in the broad conversion region. The conversion probability is periodic in fiB since the 
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precession length is dependent on fiB, which determines the phase at the outer edge of the 
isotopically neutral region. The conversion probability for a low metallicity, 25 M model 
is not shown here. Like for the 15M model conversion is effective for a field of Bq = 10 10 
G but it contains a number of complicated features due to precession in the inner part of 
the isotopically neutral region. 

Also it is found that in contrast to the metallicity the dependence on the progenitor mass 



is comparatively weak |10 



IV. STATISTICAL ANALYSIS 

In total 19 neutrinos from SN1987A have been detected: Kamiokande recorded 11 and 
1MB eight events A maximum likelihood analysis is frequently used as an effective and 
unbiased method to interpret this sparse data [P,^4[. I will be using the likelihood function 
suggested by Jegerlehner et al. |J. If the expected spectrum of detected energies is n(E), 
the total number of events N ohs and the observed energies E iy the likelihood function C is 
given by 

/ />00 \ -^obs 

£ oc exp n(E)dEj J[ n(Ei). (6) 

For a joint analysis of the signal from both detectors, the likelihood function is the product 
of the individual £'s. The expected spectrum of detected energies is given by 

/■oo 

n(E)= dE'P(E, E')rj(E')n p (E'), (7) 
Jo 

where P denotes the probability of detecting an energy E if the real energy is E', which is 
assumed to be Gaussian with an energy-dependent dispersion: 

where E s = 0.75 MeV for Kamiokande and E s = 1.35 MeV for 1MB. rj represents the 
efficiency curve for the respective detector. Analytic fits to these curves |TJ for Kamiokande 
are 



f)(E) = 0.93 - exp[-(£/9.0 MeV) 2 ' 5 ] (9) 

for E > 7.0 MeV and rj = for E < 7.0 MeV. For the 1MB detector 

r)(E) = 0.3975 — ^— - 0.02625 ( — - 0.59 (10) 
/v ; 10 MeV V 10 MeV/ v ; 

for 1.9 < E/10 MeV < 7.6, r] = 0.915 for E/10 MeV > 7.6 and 77 = for £/10 MeV < 1.9. 
Lastly 

N 

n P (E) = ^-A E + Q) F ^ E + Q)» ( n ) 

where Q is the mass difference between neutron and proton, Q = 1.29 MeV, a = 2.295 x 
10~ 44 cm 2 is the cross-section for the superallowed reaction u e (p, n)e + , F Pe is the electron 
antineutrino flux at the detector, D = 50 kpc is the distance to the supernova and N p 
is the number of target protons in the detector: N p = 1.43 x 10 32 for Kamiokande and 
N p = 4.55 x 10 32 for 1MB. 

The likelihood function is a function of the free parameters of the model Xi and the data 
Hi. The best-fit values x~i are the parameters which maximize the likelihood function 

£(x, y) = maxC(x, y) (12) 

Under the assumption that C is Gaussian, the confidence region around the best fit values 
Xi can be estimated by 

ln£(x,y)-ln£(x,y)<^ X , (13) 

where x is 2.3, 4.61 and 6.17 for the 68.3%, 90% and 95.4% confidence levels respectively. 

In the presence of SF-oscillations (P e <-> u^) the primary electron anti-neutrino spectrum 
becomes distorted due to the admixture of the spectrum, so that the flux at some distance 
from the source is related to the primary spectra F° by 

F u - e = (l-p)F° Pe +pF^ (14) 
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where the conversion probability p is a function of the neutrino energy E and the mixing pa- 
rameters. The time-integrated primary spectra are approximated by a Maxwell-Boltzmann 
distribution 

F Q V {E) = CE 2 e~ E ^ (15) 

with the parameters C and T v . Although there is little reason to believe that the actual 
spectra really follow a Maxwell-Boltzmann distribution, it allows to fit the first and second 
moments, < E > and < E 2 >, of the spectra and seems reasonable in the light of an 
exponential cooling model. 

Equipartition of energy between all (anti-) neutrino species is assumed, but different 
neutrino 'temperatures' for the v e s and v^s, T De and T u , are allowed. Due to their lower 
interaction cross-section the temperature of the v^s is higher and I define 

T v>i = aT De (16) 

with a predicted to lie between 1.5 and 2.0. The total energy E t is six times the energy 
emitted in v e s which is found to be 

E t = 36CT 4 . (17) 

As in Ref. || the detector background is expected to have a negligible influence on the 
results and is ignored here. 

V. NUMERICAL RESULTS 

Using the method described above the best-fit values for T Pe and E t were first found 
in the absence of oscillations. The Kamiokande data yielded a value for T 9e of 2.6 MeV 
with a total energy of E t = 4.8 x 10 53 erg while the 1MB data gave T Pe = 3.6 MeV and 
E t = 4.7x 10 53 erg. These results are in good agreement with the values found by Jegerlehner 
et al. [§ of T Pe = 2.5 MeV and E t = 4.9 x 10 53 erg for Kamiokande and T De = 3.7 MeV and 
E t = 5.4 x 10 53 erg for 1MB and thus independently confirm their results. 

9 



Subsequently I calculated the best-fit parameters for a variety of models. Best fit values 
for E t and T Pe together with errorbars indicating the 90% confidence levels are shown in Fig. 
3 and 4 for the Kamiokande and the 1MB data respectively. For the Z & - models the mixing 
is too small to have a significant effect and the results are virtually indistinguishable from 
the ones obtained in the absence of oscillations and not shown in Fig. 3 and 4. However the 
strong conversion in the Zq- models shifts the inferred T Pe to significantly lower values. Also 
a higher relative temperature of the muon neutrinos (i.e. a higher a) lowers the inferred 
neutrino temperature. For the Kamiokande data T Pe is lowered to 1.9 MeV and 1.5 MeV for 
a = 1.5 and 2.0 respectively. Similarly for the 1MB set the inferred value for T 9e is 2.5 MeV 
and 1.8 MeV. 

The overlap between the values deduced from the two data sets is poor and I refrain 
from a joint analysis of the data. 

VI. CONCLUSION 

In agreement with other authors it is found that in the absence of oscillations the results 
for the total energy released in neutrinos are in excellent concordance with supernova theory 
(1), whereas the temperature T Pe lies slightly below the predicted range (2). The lower end 
of the predicted range lies within the 90% confidence region. 

As expected any admixture of the hotter spectrum through SF-oscillations aggravates 
the discrepancy between the observed and predicted neutrino spectra. It is seen that a bound 
on the neutrino transition magnetic moment fi v which is deduced from the observation of 
SN neutrinos is strongly dependent on the metallicity of the progenitor. It is found that 
a low metallicity can substantially increase the sensitivity to \x and has to be taken into 
account when deducing bounds on \i v from supernova neutrinos. Therefore bounds on \x v 
derived from the neutrino signal from SN1987A may be less stringent than earlier thought 
0. For solar metallicities the data from SN1987A is consistent with a magnetic moment of 
fi u ~ 10~ 12 [iB for an y realistic magnetic field inside the supernova and mixing parameters 
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chosen from the MSW solution to the solar neutrino puzzle. However lower metallicities as 
found in the LMC make a magnetic moment as high as lO~ 12 fj, B seem less likely unless the 
magnetic field is lower than argued here. 

Thus under the assumption that the progenitor of SN1987A in fact had a metallicity 
as low as Z and a magnetic field as high as B « 10 10 G the data seems to restrict the 
magnetic moment to \i v < lCr 13 /x#. 

The aim of this article is to point out the dependence of the inferred value for the neutrino 
magnetic moment on the metallicity of the progenitor thus highlighting the uncertainties 
involved in interpreting the data. I should emphasize again that the poor statistics of the 
data as well as the many unknowns of the source do not allow a reliable determination of 
neutrino properties such as the transition magnetic moment. However there is hope that 
with the advent of a new generation of neutrino telescopes as Superkamiokande, the Large 
Volume Detector (LVD), the Sudbury Neutrino Observatory (SNO) and others a future 
supernova explosions in our neighbourhood will provide us with richer data allowing better 
estimates of neutrino properties and variables of the supernova explosion. 
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FIGURES 

FIG. 1. Magnitude of the effective potential AV as a function of radius for several progenitor 
masses and metallicities. AV is measured in eV and the radius is measured in units of solar radii, 
(a) 15 M , Z , (b) 15 M , Z , (c) 25 M , Z , (d) 25 M , Z 

FIG. 2. Conversion probabilities as function of neutrino energy E u , solid line: Z , 15 M , 
B = 10 10 G, dotted line: Z , 25 M , B = 10 10 G, dashed line: Z , 15 M , B = 10 10 G, 
dot-dashed line: Z , 25 M , B = 10 9 G 

FIG. 3. Best fit values for E t and T Ve for the Kamiokande data. The errorbars correspond to 
90% confidence level. 

FIG. 4. The same as Fig. 3 but for the 1MB data. 
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